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PULSE OXIMETER WITH ALTERNATE HEART-RATE 

DETERMINATION 

CROSS-REFERENCES TO RELATED APPLICATIONS 
5 [00011 NOT APPLICABLE 

STATEMENT AS TO RIGHTS TO INVENTIONS MADE UNDER 
FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT 
[0002] NOT APPLICABLE 

10 

REFERENCE TO A "SEQUENCE LISTING," A TABLE, OR A COMPUTER 
PROGRAM LISTING APPENDIX SUBMITTED ON A COMPACT DISK 
[0003] NOT APPLICABLE 

1 5 BACKGROUND OF THE INVENTION 

[0004] The present invention relates to oximeters, and in particular to determining a pulse 
rate by multiple mechanisms in a detected waveform from a pulse oximeter. 

[0005] Pulse oximetry is typically used to measure various blood chemistry characteristics 
including, but not limited to, the blood-oxygen saturation of hemoglobin in arterial blood, the 

20 volume of individual blood pulsations supplying the tissue, and the rate of blood pulsations 
corresponding to each heartbeat of a patient. Measurement of these characteristics has been 
accomplished by use of a non-invasive sensor which scatters light through a portion of the 
patient's tissue where blood perfuses the tissue, and photoelectrically senses the absorption of 
light in such tissue. The amount of light absorbed at various wavelengths is then used to 

25 calculate the amount of blood constituent being measured. 

[0006] The light scattered through the tissue is selected to be of one or more wavelengths 
that are absorbed by the blood in an amount representative of the amount of the blood 
constituent present in the blood. The amount of transmitted light scattered through the tissue 
will vary in accordance with the changing amount of blood constituent in the tissue and the 
30 related light absorption. For measuring blood oxygen level, such sensors have typically been 



provided with a light source that is adapted to generate Hght of at least two different 
wavelengths, and with photodetectors sensitive to both of those wavelengths, in accordance 
with known techniques for measuring blood oxygen saturation. 

[0007] Known non-invasive sensors include devices that are secured to a portion of the 
5 body, such as a finger, an ear or the scalp. In animals and humans, the tissue of these body 
portions is perfused with blood and the tissue surface is readily accessible to the sensor. 

[0008] Patents No. 6,083,172, No. 5,853,364 and No. 6,41 1,833 show multiple methods of 
calculating a pulse rate in a pulse oximeter, with a "best rate" module which arbitrates 
between the pulse rate calculations to select a best rate based on confidence levels associated 
10 with each. The confidence levels are calculated using various metrics to determine the 

reliability of the different pulse rate calculations. Also, Patent No. 5,524,631 shows a fetal 
heart rate monitor that uses multiple parallel filter paths to identify the fetal heart rate, and 
uses a figure of merit operation to weight the different heart rate estimates. 

[0009] N-100 . The N-lOO technology, dating to around 1985, accepted or rejected pulses 
15 based on pulse history of the size of pulses, pulse shape, expected time to occur (frequency) 
and ratio of R/IR. 

[0010] In particular, the N-100 found pulses by looking for a signal maximum, followed by 
a point of maximum negative slope, then a minimum. The processing was done in a state 
machine referred to as "munch." Each maximum was not qualified until the signal passed 

20 below a noise threshold, referred to as a noise gate. This acted as an adaptive filter since the 
noise gate level was set by feedback fi-om a subsequent processing step to adapt to different 
expected signal amplitudes. The pulses are then accepted or rejected in a "Level3" process 
which was a filter which adapts to changing signals by comparing the amplitude, period and 
ratio-of-ratios (ratio of Red to IR, with Red and IR being expressed as a ratio of AC to DC) of 

25 a new pulse to the mean of values in a history buffer, then determining if the difference is 

within a confidence level. If the new pulse was accepted, the history buffer was updated with 
the values for the new pulse. The leveB process acted as an adaptive bandpass filter with 
center- frequency and bandwidth (confidence limits) being adapted by feedback from the 
output of the filter. 

30 [0011] N-2Q0 . The N-200 improved on the N-100 since it could be synchronized with an 
ECG, and included ECG filtering. The N-200 also added interpolation to compensate for 
baseline shift between the time of measuring the pulse maximum and minimum. The N-200 
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included other filtering features as well, such as a "boxcar" filter which computed the mean 
of a varying number of signal samples. 

[0012] The N-200, after various filtering and scaling steps, applies the digitized signals to a 
"boxcar" filter, which computes the mean of N samples, where N is set by feedback fi-om a 
5 subsequent processing step according to the filtered heart rate. New samples are averaged 
into the boxcar filter, while the oldest samples are dropped. The boxcar length (N) is used to 
set three parameters: a pulse threshold, absolute minimum pulse and small pulse. An 
ensemble-averaging (a.k.a "slider") filter then produces a weighted average of the new 
samples and the previous ensemble-averaged sample firom one pulse-period earlier. The 
10 samples are then passed to a "munch" state machine and a noise gate, like the N-100. An 

interpolation feature is added to the N-100 process, to compensate for changes in the baseline 
level. Since the minimum and maximum occur at different times, a changing baseline may 
increase or decrease the minimum and not the maximum, or vice-versa, 

[0013] "Ensemble averaging" is an integral part of C-Lock, which is NELLCOR's 
1 5 trademark for the process of averaging samples firom multiple pulses together to form a 
composite pulse. This process is also known as "cardiac-gated averaging." It requires a 
"trigger" event to mark the start of each pulse. 

BRIEF SUMMARY OF THE INVENTION 
20 [0014] The present invention is a pulse oximeter which determines multiple heart rates, and 
selects between them based on the metrics of only one of the heart rate calculations. A 
primary heart rate calculation method is selected, and is used unless its metrics indicate 
questionable accuracy, in which case an alternative rate calculation is available and is used 
instead. 

25 [0015] In one embodiment, the primary heart rate calculation method does not use an 
ensemble averaged waveform, while the alternative heart rate calculation does use an 
ensemble averaged waveform. The altemative heart rate calculation is used if the primary 
calculation has disqualified its most recently detected pulse. 

30 BRIEF DESCRIPTION OF THE DRAWINGS 

[0016] Fig. 1 is a block diagram of an oximetry system incorporating an embodiment of the 
invention. 
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[0017] Fig. 2 is a diagram of the software processing blocks of an oximeter including an 
embodiment of the present invention. 

[0018] Fig. 3 is a context diagram of the pulse rate calculation subsystem. 



5 DETAILED DESCRIPTION OF THE INVENTION 

[0019] Fig. 1 illustrates an embodiment of an oximetry system incorporating the present 
invention. A sensor 10 includes red and infrared LEDs and a photodetector. These are 
connected by a cable 12 to a board 14. LED drive current is provided by an LED drive 
interface 16. The received photocurrent from the sensor is provided to an I-V interface 18. 

1 0 The IR and red voltages are then provided to a sigma-delta interface 20 incorporating the 

present invention. The output of sigma-delta interface 20 is provided to a microcontroller 22 
which includes a 10-bit A/D converter. Controller 22 includes flash memory for a program, 
and EEPROM memory for data. The processor also includes a controller chip 24 connected 
to a flash memory 26. Finally, a clock 28 is used and an interface 30 to a digital calibration 

15 in the sensor 10 is provided. A separate host 32 receives the processed information, as well 
as receiving an analog signal on a line 34 for providing an analog display. 

[0020] Design Summary The design of the present invention is intended to deal with 
unwanted noise. Signal metrics are measured and used to determine filter weighting. Signal 
metrics are things that indicate if a pulse is likely a plethysmograph or noise, such as 
20 frequency (is it in the range of a human heart rate), shape (is it shaped like a heart pulse), rise 
time, etc. A similar technique was used in the Nellcor N200, described in the background of 
this application. The new design adds a number of different features and variations, such as 
the use of two ensemble averagers as claimed in the present invention. 

[0021] Details of the architecture are shown in the diagram of Fig. 2. This design 
25 calculates both the oxygen saturation, and the pulse rate, which are described separately 
below. 

[0022] L Oxvgen Saturation Calculation. 

[0023} A. Signal Conditioning - The digitized red and IR signals are received and are 
conditioned in this block by (1) taking the 1st derivative to get rid of baseline shift, (2) low 
30 pass filtering with fixed coefficients, and (3) dividing by a DC value to preserve the ratio. 
The function of the Signal Conditioning subsystem is to emphasize the higher frequencies 
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that occur in the human plethysmograph and to attenuate low frequencies in which motion 
artifact is usually concentrated. The Signal Conditioning subsystem selects its filter 
coefficients (wide or narrow band) based on hardware characteristics identified during 
initialization. 

5 [0024] Inputs - digitized red and IR signals 

[0025] Outputs - Pre-processed red and IR signals 

[0026] B. Pulse Identification and Qualification - The low pass filtered and digitized red 
and IR signals are provided to this block to identify pulses, and qualify them as likely arterial 
pulses. This is done using a pre- trained neural net, and is primarily done on the IR signal. 
1 0 The pulse is identified by examining its amplitude, shape and frequency, just as was done in 
the Nellcor N-100. An input to this block is the average pulse period from block D. This 
function is similar to the N-100, which changed the upfi-ont qualification using the pulse rate. 
The output indicates the degree of arrhythmia and individual pulse quality. 

[0027] Inputs - (1) Pre-processed red and IR signals, (2) Ave. pulse period, (3) 

15 Lowpass Waveforms from the low pass filter. 

[0028] Outputs - (1) Degree of arrhythmia, (2) pulse amplitude variations, (3) 

individual pulse quality, (4) Pulse beep notification, (5) qualified pulse periods and age. 

[0029] C. Compute Signal Quality Metrics - This block determines the pulse shape 
(derivative skew), period variability, pulse amplitude and variability. Ratio of Ratios 
20 variability, and frequency content relative to pulse rate. 

[0030] Inputs - (1) raw digitized red and IR signals, (2) degree of arrhythmia, 

individual pulse quality, pulse amplitude variation (3) pre-processed red and IR signals, 
(4) average pulse period. 

[0031] Outputs - (1) Lowpass and ensemble averaging filter weights, (2) metrics for 

25 sensor off detector, (3) Normalized Pre-processed waveforms, (4) percent modulation. 

[0032] D. Average Pulse Periods. This block calculates the average pulse period from the 
pulses received. 

[0033] Inputs - Qualified pulse periods and age. 

[0034] Outputs - Average pulse period. 
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[0035] El . Lowpass Filter and Ensemble Averaging - Block El low pass filters and 
ensemble averages the signal conditioned by block A, and normalized by block C, for the 
pulse rate identification. The weights for the low pass filter are determined by the Signal 
Metrics block C. The signal is also ensemble averaged (this attenuates frequencies other than 
5 those of interest near the pulse rate and its harmonics), with the ensemble averaging filter 
weights also determined by Signal Metrics block C. Less weight is assigned if the signal is 
flagged as degraded. More weight is assigned if the signal is flagged as arrhythmic because 
ensemble-averaging is not appropriate during arrhythmia. Red and LR are processed 
separately, but with the same filtering weights. The filtering is delayed approximately one 
10 second to allow the signal metrics to be calculated first. 

[0036] The filters use continuously variable weights. If samples are not to be ensemble- 
averaged, then the weighting for the previous filtered samples is set to zero in the weighted 
average, and the new samples are still processed through the code. This block tracks the age 
of the signal - the accumulated amount of fihering (sum of response times and delays in 
1 5 processing). Too old a result will be flagged (if good pulses haven't been detected for 
awhile). 

[0037] Inputs - (1) normalized pre-processed red and IR signals, (2) average pulse 

period, (3) low pass filter weights and ensemble averaging filter weights, (4) ECG triggers, if 
available, (5) IR fundamental, for zero-crossing triggers. 

20 [0038] Outputs - (1) filtered red and IR signals, (2) age. 

[0039] F. Estimate Filtered Waveform Correlation and Calculate Averaging Weight - this 
uses a noise metric similar to that used in the NlOO and N200 described above, and doesn't 
use feedback. The variable weighting for the filter is controlled by the ratio-of-ratios 
variance. The effect of this variable-weight filtering is that the ratio-of-ratios changes slowly 
25 as artifact increases and changes quickly as artifact decreases. The subsystem has two 

response modes. Filtering in the Fast Mode targets an age metric of 3 seconds. The target age 
is 5 seconds in Normal Mode. In Fast Mode, the minimum weighting of the current value is 
clipped at a higher level. In other words, a low weight is assigned to the newest ratio-of- 
ratios calculation if there is noise present, and a high weight if no noise is present. 

30 [0040] Inputs - (1) filtered red and IR signals and age, (2) calibration coefficients, (3) 

response mode (user speed settings). 
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[0041] 



Outputs - averaging weight for ratio-of-ratios calculation. 



[00421 H. Calculate Saturation - Saturation is calculated using an algorithm with the 
calibration coefficients and averaged ratio of ratios. 



[0043] 



Inputs - (1) Averaged Ratio-of-Ratios, (2) calibration coefficients. 
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[0044J 



Outputs - Saturation. 



[0045] IL Pulse Rate Calculation. 

[0046] E2. Lowpass Filter and Ensemble Averaging - Block E2 low pass filters and 
ensemble averages the signal conditioned by block A, for the pulse rate identification. The 
weights for the low pass filter are determined by the Signal Metrics block C. The signal is 

10 also ensemble averaged (this attenuates frequencies other than those of interest near the pulse 
rate and its harmonics), with the ensemble averaging filter weights also determined by Signal 
Metrics block C. Less weight is assigned if the signal is flagged as degraded. More weight is 
assigned if the signal is flagged as arrhythmic since filtering is not appropriate during 
arrhythmia. Red and JR are processed separately. The process of this block is delayed 

15 approximately one second to allow the signal metrics to be calculated first. 

[0047] The filters use continuously variable weights. If samples are not to be ensemble- 
averaged, then the weighting for the previous filtered samples is set to zero in the weighted 
average, and the new samples are still processed through the code. This block tracks the age 
of the signal - the accumulated amount of fiUering (sum of response times and delays in 
20 processing). Too old a result will be flagged (if good pulses haven't been detected for 
awhile), 

[0048] Inputs - (1) pre-processed red and IR signals, (2) average pulse period, (3) 

Lowpass filter weights and ensemble averaging filter weights, (4) ECG triggers, if available, 
(5) IR fundamental, for zero-crossing triggers. 

25 [0049] Outputs - (1) filtered red and IR signals, (2) age. 

[0050] I. Filtered Pulse Identification and Qualification - This block identifies and 
qualifies pulse periods from the filtered waveforms, and its results are used only when a pulse 
is disqualified by block B. 

[0051] Inputs - (1) filtered red and IR signals and age, (2) average pulse period, (3) 

30 hardware ID or noise floor, (4) kind of sensor. 
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[0052] Outputs - qualified pulse periods and age. 

[0053] J. Average Pulse Periods and Calculate Pulse Rate - This block calculates the pulse 
rate and average pulse period. 

[0054] Inputs - Qualified pulse periods and age 

5 [0055] Outputs - (1) average pulse period, (2) pulse rate. 

[0056] III. Venous Pulsation 

[0057] K. Detect Venous Pulsation - Block K receives as inputs the pre-processed red and 



IR signal and age from Block A, and pulse rate and provides an indication of venous 
pulsation as an output. This subsystem produces an IR fundamental waveform in the time 
10 domain using a single-tooth comb filter which is output to the Ensemble Averaging filters. 

[0058] Inputs - (1) filtered red and IR signals and age, (2) pulse rate. 

[0059] Outputs - Venous Pulsation Indication, IR fundamental 

[0060] IV. Sensor Off 

[0061] L. Detect Sensor-Off and Loss of Pulse Amplitude - The Pulse Lost and Sensor Off 
15 Detection subsystem uses a pre-trained neural net to determine whether the sensor is off the 
patient. The inputs to the neural net are metrics that quantify several aspects of the behavior 
of the IR and Red values over the last several seconds. Samples are ignored by many of the 
oximetry algorithm's subsystems while the Signal State is not either Pulse Present or Sensor 
Maybe Off. The values of the Signal State variable are: "Pulse Present, Disconnect, Pulse 
20 Lost, Sensor Maybe Off, and Sensor Off." 

[0062] Inputs - (1) metrics, (2) front-end servo settings and ID 

[0063] Outputs - Signal state including sensor-off indication 

Pulse Rate Calculation subsystem 

[0064] The subsystem averages qualified pulse periods from the Pulse Identification and 
25 Qualification subsystem. It outputs the average period and the corresponding pulse rate. 

[0065] The oximetry algorithm contains two instances of this subsystem. The first instance 
receives input from the Pulse Identification and Qualification instance whose input waveform 
have been processed by the Signal Conditioning subsystem, then lowpass filtered, but not 
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ensemble averaged, by the Ensemble Averaging subsystem. The second instance of the Pulse 
Rate Calculation subsystem receives input from two instances of the Pulse Identification and 
Qualification subsystem, the one described above and a second instance that receives input 
that has been ensemble averaged. 



5 Selection of Pulse Period Source 

[0066] One instance of the subsystem receives qualified pulse periods from two sources. 
The subsystem selects which of these two sources to use for its pulse rate calculation based 
solely on analysis of only one source, the "primary" source. The oximetry algorithm 
designates the Pulse Identification and Qualification instance that does NOT receive 

1 0 ensemble-averaged waveforms as the primary source, and designates the other Pulse 

Identification and Qualification instance as the "altemate" source of qualified pulse periods. 
Qualified pulse periods from the altemate source are only used if the most recent pulse from 
the primary source was rejected. When a qualified pulse period is received from the primary 
source, it is always used to update the pulse-rate calculation, and will prevent qualified pulse 

15 periods from the altemate source from being used until the primary source once again rejects 
a pulse period. 

Calculation of Average Pulse Period and Pulse-Rate Estimate 

[0067] When the subsystem uses a Qualified_Pulse_Period, it updates its average pulse 
period, Avg_Period, using a pulse-based, variable-weight IIR filter, then computes its Rate 
20 output from Avg Period. The steps for this filtering operation are: 

1. K=^(60/At)/ Qualified_Pulse_Period 

2. k= ConsecutivejQualified I max(|r, - |r,.i - r,_2|, \r,.2 - ^-sIj 1-0) 

3. x = hound{m\n{Avg_Periodt -i, Qualified_Pulse_Period), Va seconds, 2 seconds) / 7 seconds 

4. If Rate_Age > 10 seconds, x = min(x * Rate_Age / 10 seconds, 0.3) 
25 5. k= max(l / Total jQualified, mm{k, x)) 

6. If Avg_Periodt .\ <> 0 Avg_Periodt = Avg_Periodt a * {Qualified_Pulse_Period I 
Avg_Periodt 

7. If Avg_Periodt -i = 0 Avg_Periodt = Qualified_Pulse_Period 

8. Rate = (60 / AO / Avg_Periodt 

30 9. Rate_Age = Rate_Age + A: * {Qualified_Penod_Age - Rate_Age) 
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where: 

r, is the pulse rate corresponding to Qualified_Pulse_Period, in BPM 
the ,.| subscript denotes the previous qualified pulse. 
A/ is the oximetry algorithm's sample interval in seconds 
5 60 / At is the number of samples per minute 

jc is a filter weight that targets a 7-second response time for typical adult pulse rates. 
k is the final filter weight, based on both x and the differences between consecutive values of 
ri. During the first few pulses, k is increased to at least 1 / Total jQualified so that the initial 
qualified pulses will be weighted equally. 

1 0 [0068] ConsecutivejQualified is the number of consecutive qualified pulses, and 

Total jQualified is the total number of pulses qualified since the subsystem was reinitialized. 
Both ConsecutivejQualified and Total JQualified are incremented each time a 
Qualified_PulseJPeriod is used, before k is calculated. ConsecutivejQualified is set to zero 
when a pulse is rejected by the pulse-period source currently in use. 

1 5 [0069] The update formula for Avg^Periodt^ in step 6 above, is a geometric average of 

Avg_Periodt and Qualifiied_PulseJPeriod. Geometric averaging helps to keep the subsystem 
responsive to large pulses-to-pulse period variations, and large, sustained changes in pulse 
rate. 

[0070] Once Rate is initialized to a non-zero value, Rate_Age is incremented every sample, 
20 whether or not Rate is updated. 

Context Diagram 

[0071] Fig. 3 is a context diagram of the pulse rate calculation subsystem. The subsystem 
updates its Avg_Period and Rate outputs fi'om Qualified JPulse_Periods, It uses 
Qualified PulseJPeriods from the Alternative J^eriodJSource only if it last received a 
25 Notify _Pulse_Rejected fi"om the primary source. It updates its Rate Age output based on 
Qualified_Period_Age, When Rate is updated, the subsystem sets its Pulse_Rate Updated 
flag. The Reinitialize input tells the subsystem to reinitialize itself. Increment JR.ate_Age 
notifies the subsystem to increment its Ratejlge every sample once Rate is initialized. 
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